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The Fe center in ferroelectric PbTi03 together with an oxygen vacancy forms a charged defect 
associate, oriented along the crystallographic c-axis. Its microscopic structure has been analyzed 
in detail comparing results from a semi-empirical Newman superposition model analysis based on 
finestructure data and from calculations using density functional theory. 

Both methods give evidence for a substitution of Fe 3+ for Ti 4+ as an acceptor center. The position 
of the iron ion in the ferroelectric phase is found to be similar to the B-site in the paraelectric phase. 
Partial charge compensation is locally provided by a directly coordinated oxygen vacancy. 

Using high-resolution synchrotron powder diffraction, it was verified that lead titanate remains 
tetragonal down to 12 K, exhibiting a c/o-ratio of 1.0721. 

Keywords: lead titanate ceramics, ferroelectrics, iron functional center, oxygen vacancy, charge compensa- 
tion, EPR, Newman superposition model, DFT calculations 
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I. INTRODUCTION 

Piezoelectric lead titanate oxide (PbTiC>3, PT) is of 
considerable scientific and technical interest because of 
its extraordinary electromechanical properties leading to 
the development of novel devices, such as nonvolatile 
memories, detectors, sensors or actuatorsi 2 ^. A variety 
of ferroelectric properties can be controlled by replacing 
the lead or titanium cations with rare-earth or transition- 
metal ions. In general, the effect of aliovalent dopants on 
PbTiC>3, such as Fe 3+ , leads to the creation of oxygen 
vacancies (V") (The Kroger and Vink notation is used 
to designate the charge state of the defect with respect to 
the neutral lattice.) for charge compensation, which are 
expected to have a major impact on properties and per- 
formance of ferroelectric compounds. However, it is still 
not certain at which coordination sphere the charge com- 
pensation takes place, and because the oxygen vacancies 
are known to be the dominant charge carriers in this class 
of compounds this is an issue of technological and scien- 
tific importance 4 . Beside charge compensation provided 
by V" , intrinsic doubly negatively charged Vp h centers 
have been proposed as additional charge compensating 
mechanism^, these vacancies, however, being rather im- 
mobile and almost of no importance for charge transport. 

When dealing with dopants on a sub-percentage level, 
a sensitive test of the local environment around the func- 
tional center can only be provided by electron param- 
agnetic resonance (EPR). In case of paramagnetic Fe 3+ 
dopants, the sextet spin ground state is most influenced 
by the ligand-field originating from the nearest-neighbor 
2 ~ ions. The resulting finestructure (FS) interaction is 
a measure of the local symmetry. Various structure mod- 
els have been proposed such as Fe' Ti — V" defect asso- 
ciates, not-coordinated 'free' Fe^j centers, as well as off- 
center shifted iron centers. Previous EPR studies on crys- 
talline samples reported FS parameters spanning a large 



range 6 . Structural information about the dopant site in 
PbTiC>3 and related compounds can be based on mod- 
elling the magnitude of the FS parameter, and therefore 
its accurate determination is an important issue. Because 
most samples are provided in polycrystalline form, it was 
important to demonstrate that FS parameters can reli- 
ably be determined even for polycrystalline compounds 
by invoking high-frequency-EPRt 

Structure modelling can be performed using different 
levels of sophistication. For instance, the FS parameter 
can be analyzed in terms of the semi-empirical Newman 
superposition model (NSM) S . Within this model, values 
for local distortion, i.e., atomic displacements from the 
ideal crystal structure in the vicinity of the Fe 3+ cen- 
ter, are derived by comparing calculated and observed 
FS parameters. A detailed study was performed already 
25 years ago by Siegel and Mulle»2*iS, in which the dis- 
tortion of the cubic high-temperature unit cell in the 
ferroelectric phase was taken as decisive parameter con- 
trolling the FS splitting at the iron site. Over the past 
decades, the NSM parametrization has been refined by 
applying the method to many examples of paramagnetic 
centers in perovskite-type crystals. Examples include 
iron-doped PbTi0 3 , SrTi0 3 and BaTiOj^iS, LiTaOgii, 
as well as chromium-modified PbTiOg^ and manganese 
in BaTiOgi 3 .. 

Second, an alternate approach offering additional mi- 
croscopic information, is provided by calculating the 
equilibrium structure using density-functional theory 
(DFT). In recent years numerous calculations have 
been performed on ferroelectric materials in order to 
explain the microscopic mechanisms of spontaneous 
polarizationiiii^ and phase diagramsi§*i£, or to deter- 
mine the full piezoelectric stress tensor and dynamical 
charges^. In particular, the role of defects and defect 
dipoles has been highlighted. Oxygen vacancies were 
shown to pin at 180°-domain walls, confirming the ten- 
dency of these defects to migrate to the walls 21 . It was 
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can be written as 



Figure 1: Schematic representation of the tetragonal unit cell 
of PbTiOa with designation of the different oxygen sites. The 
experimental atomic distances for the Ti-site to the surround- 
ing oxygen octaeder are given in pm. 



also derived that lead vacancies have stable charges rang- 
ing from —2 to —4 and thus being an effective accep- 
tor in lead titanate22. The dipole moment V" — Vp h of 
di-vacancies was calculated, demonstrating that the in- 
creased value may be an important source of local po- 
larization and electric fields^ 3 -. Another defect dipole 
originates from interstitial hydrogen impurities in lead 
titanate, which were found to bind to oxygen and to act 
as shallow donor impurities. The H-0 dipole increases 
the polarization and the barrier for reversing the defect 
dipole, which can give rise to imprints^*. 

In this work, we exploit the recently accurately deter- 
mined FS parameter 7 - as basis for structure modelling by 
performing an NSM analysis. Second, we compare the 
results obtained from this semi-empirical approach with 
results from DFT calculations, thus providing a sound 
basis for the postulated structural relaxation at the im- 
purity site. As a result, we present a refined structural 
model of the iron functional center in lead titanate. 



II. THEORETICAL METHODS 



A. Newman superposition model 



With its five unpaired electrons, the free Fe 3+ ion pos- 
sesses a half-filled 3d shell and can be described as orbital 
singlet. Its ground state configuration is 6 £5/2 (S = §) 
and the six-fold spin degeneracy can be lifted by the FS 
interaction and an external magnetic field. Neglecting 
hyperfine interaction, because the magnetic active iso- 
tope 57 Fe is present in only 2.2 % natural abundance, an 
approximate spin-Hamiltonian for this high-spin system 



H = (3 e B • g • s 
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Here, (3 e denotes the Bohr magneton, Bo is 
ternal magnetic field, g the electron g-matrix, 
the FS Hamiltonian coefficients and 0\ are the extended 
Stevens spin operators2^S&. The first term represents the 
electronic Zeeman interaction and the second term is the 
effective FS Hamiltonian, describing the interaction of 
the crystal field with the paramagnetic ion. The rank k 
in the Hamiltonian must be even, and k is restricted by 
k < 2S and q < k, allowing terms up to k = 4 for S = |. 
However, parameters b q k only up to second rank (fc = 2) 
will be used, because terms of fourth rank were shown to 
be at least two orders of magnitude smaller than those of 
second rank&(c-e). Parameters b q k are related to the stan- 
dard FS parameters B q k by means of the scaling factor as 
B k = fkb q k - For k — 2 one finds f% = \. The coefficients 
bl are related to the conventional spectroscopic FS pa- 
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The FS term in lifts the degeneracy of eigen states 
even in the absence of an external magnetic field. There 
is always a choice of spin-operator coordinate system with 
respect to the crystal axes for which the ratio b\fb\ lies in 
the interval [0, 1], analogous to the condition E/D < 1/3. 

In case the symmetry of the paramagnetic site is axial, 
the FS coefficients in its eigen frame reduce to b® 0> 
b\ = 0. Furthermore, the g-matrix for S-state ions usu- 
ally has very small anisotropy and can be treated as 
isotropic. The resulting scalar g value is expected to be 
close to that of the free electron g e . For Fe 3+ we therefore 
set <7i so = 2.002 and then the spin-Hamiltonian used for 
numerical spectrum simulation of the experimental data 
reduces to 



n = A, 9iso B • s + bl 



s 2 z - \s{s + i) 



(2) 



The Newman superposition mode? allows for an anal- 
ysis of positions of neighboring ions (ligands) around the 
Fe 3+ center, given the single-ligand contributions b k (Ri) 
are known. The essential assumption of the NSM is that 
the spin-Hamiltonian parameters for a paramagnetic ion 
can be constructed by a superposition from individual 
contributions of separate neighboring ligands. The con- 
tribution of next nearest neighbor ions as well as interac- 
tion between the ligands are ignored. The NSM expres- 
sion for the zero-field splitting parameters can be formu- 
lated as follows 



bl = J2h(Ri)Kl( 



(3) 



Here, i denotes every ligand, b k is the contribution 
from each single ligand, Ri is the distance between the 



B Density functional theory (DFT) modelling 

ith ligand and the paramagnetic ion, Oi and <f>i are the 
polar and axial angles between i?j and the symmetry axis 
of the paramagnetic center, and K%(9i, <j>i) are spherical 
harmonic functions of rank k of the polar angles, listed 
in£(f). For the axial second-rank parameter one has 

K° 2 =\(i cos 2 9-1) (4) 

For the short-range distance dependence of the single- 
ligand contribution, the following empirical power law 
was established 

b k (Ri) = h(R ) (J^j k (5) 

The critical exponent parameter tk is specific to a par- 
ticular ion-ligand system. For a given ligand, the intrin- 
sic parameters bk(Ri) are determined by the nature of the 
ligand and the covalency of the bond, which is assumed 
to depend exclusively on the bond lengths Ri. Because 
bk(Ri) in this model depends only on the ligand and its 
distance, and not on other properties of the host crystal, 
this parameter can be obtained from data of the same 
ion - ligand complex measured in other host crystals. 

If the molecular coordinate system coincides with the 
principal axes system of the FS tensor, all off-diagonal 
elements vanish. For the Fe 3+ ion at a position of tetrag- 
onal site symmetry (4mm) , the NSM can be used in its 
truncated form& 




The strategy followed is to calculate the FS parameter 
by using the NSM for different conceivable structural ar- 
rangements, and to compare the calculated values with 
the experimental data. The experimental value for 6° 
has been determined reliably 7 -, and thus the interception 
points between the calculated curve and the experimen- 
tal value is used for a prediction of actually realized local 
structure. Previously, the extremum of the calculated 6§ 
dependence was used to obtain an approximate ion posi- 
tion, although a correlation between the ZFS extremum 
and the ion equilibrium position cannot be derived from 
first principles iSiiS. 

B. Density functional theory (DFT) modelling 

The electronic structure of Fe:PbTiC>3 was calculated 
using two different DFT methods. First, the local relax- 
ation of the structure around the defects, Fe-ri, and Vo, 
is studied by the Vienna Ab initio Simulation Package 
(VASP) 27 , as applied recently to study Pb-0 vacancies 
in PbTiOgS 3 .. However, instead of using ultrasoft pseu- 
dopotentials, Projector Augmented Waves (PAW) have 



3 

been applied 2 ^ 2 ^. The identical basis set has been used, 
namely Pb(5d,6s,6p), Ti(3s,3p,3d,4s,4p), and 0(2s,2p), 
as well as the LDA exchange potential 3 ^, and a 4 x 4 x 4 
k-mesh. Second, the electron density distribution p(r) 
has been computed by the NFP (new full potential) pro- 
gram package2ii22i2 3 ., which is a variant of the LMTO 
(linear muffin-tin orbitals) 34 procedure using a minimal 
basis set. 
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Figure 2: Observed and calculated diffraction profiles for 
PbTiC>3, obtained from high-resolution synchrotron powder 
diffraction at 12 K (top) together with their difference curve 
(bottom); A = 0.0699988 nm. 



III. RESULTS 
A. Low-temperature X-ray data 

The low temperature crystal structure of lead titanate 
has been the subject of some controversy over the past 50 
years. Several X-ray and optical studies as well as dielec- 
tric measurements of PT samples indicated small anoma- 
lies in the temperature regions around 173 K and 123 K. 
These anomalies were detected as changes in the neg- 
ative volume expansion coefficient in combination with 
superlattice reflections in powder samples 3 ^ below 173 K 
and as changes in lattice parameters and birefringence 
in single crystals around 183 K. They were interpreted 
as a possible second order phase transition to a different 
tetragonal, later corrected into an antiferroelectric or- 
thorhombic phased. Moreover, this earlier claimed low- 
temperature tetragonal phase was assumed to transform 
into a further tetragonal phase around 113 K 3 ^, which 
does not fit into the picture of having an orthorhombic 
phase present at this temperature range. However, nei- 
ther in neutron powder profiles^, nor in X-ray dilatomet- 
ric and optical measurements using PT single crystals, 
additional evidence was seen for a structural transition at 
low temperatures^. First-principle studies 4 " 1 - show that 
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Table I: Lattice parameters (a,b = 389.0673(11) pm, c = 416.7553(30) pm), cell volume (0.0630856(5) nm 3 ) and relative atom 
coordinates (upper rows) for PbTiC>3 at 12 K, obtained from high-resolution synchrotron powder diffraction. U% is the mean 
square of atomic displacement. Also listed are the predicted values from DFT calculations with lattice parameters a = b — 386.4 
pm, c = 404.2 pm (lower rows). Pb was kept fixed at the origin for both, Rietveld refinement and DFT calculations. 



atom 




X 


Y 


Z 


Ui ■ 100 


site symmetry 


Pb 


exp. 
DFT 


0.000000 
0.000000 


0.000000 
0.000000 


0.000000 
0.000000 


0.820(27) 


4MM(001) 


Ti 


exp. 
DFT 


0.500000 
0.500000 


0.500000 
0.500000 


0.5415(14) 
0.5346 


0.66(10) 


4MM(001) 


01 


exp. 
DFT 


0.500000 
0.500000 


0.500000 
0.500000 


0.1239(19) 
0.0904 


0.38(19) 


4MM(001) 


02 / 03 


exp. 
DFT 


0.000000 
0.000000 


0.500000 
0.500000 


0.6288(13) 
0.6032 


0.38(19) 


MM2(001) 



the thermodynamically stable phase at low temperatures 
is the tetragonal one, as all the unit-cell preserving distor- 
tions at low temperature have positive elastic constants 
and as there are no other mechanical instabilities present 
that could cause a transition to a lower symmetry group. 

To identify the low-temperature structure and to ob- 
tain reference values for the interpretation of EPR data 
recorded at low temperatures, a high-resolution syn- 
chrotron powder diffraction experiment was carried out 
in reflection geometry at B2, Hasylab in Hamburg, Ger- 
many (figure |2J). The measurement was performed at 12 
K using a He closed-cycle cryostat, equilibrating the tem- 
perature in a low-pressure helium atmosphere. The pow- 
dered sample was glued onto a Si (711) low-background 
wafer and measured using an incident beam of wave- 
length 0.0699988 nm in combination with an analyzing 
crystal and a scintillation counter—. 

For Rietveld refinement, the general structure anal- 
ysis system (GSAS)^ 3 - was used. As the peak profiles 
of the ferroelectric material are quite complex, the pro- 
files were fitted using the incorporated generalized model 
for anisotropic peak broadening^ 4 .. The structural model 
could be refined with the symmetry of P4mm, keeping 
the lead positions fixed at the origin. No supcrlattice re- 
flections were detected. The c/a - ratio is 1.0721 at this 
temperature. Lattice parameters and atom coordinates 
obtained are given in table III 61 

The detected low-temperature PAmm symmetry is 
consistent with the interpretation of the EPR spectra, for 
which a FS tensor of axial symmetry was used, thus in- 
dicating that no orthorhombic phase is present. There- 
fore it can be concluded that PbTiC>3 remains tetragonal 



down to 12 K, as already assumed earlier—. 



B. Superposition-model analysis 

The intrinsic NSM parameters ti and 62(^0) have not 
yet been determined for iron centers in lead titanate. 
However, they may be adopted from similar single crys- 
tals having equivalent Fe 3+ — O 2- bonds in octahedral co- 
ordination. For Fe 3+ in MgO at the central site with oc- 
tahedral oxygen coordination, the following set of NSM 
parameters was reported: 62 = —12.3514 GHz, ti = 8, 
i?o = 210.1 pm-. The ionic positions of the nearest oxy- 
gens were taken from the X-ray data at 12 K (cf. section 
IIII Ap . Because the site symmetry defines the direction of 
FS-tensor axes, the principal FS axis of the Fe 3+ center 
was chosen along [0 1]. 

The common reference structure for undistorted cubic 
PbTiOa is the perovskite structure, in which Ti 4+ ions 
are octahedrally coordinated by oxygen ions and are po- 
sitioned at the center of the unit cell. The Pb 2+ ions are 
twelvefold coordinated by oxygens and are located at the 
corners of the cube, whereas the O 2- ions are centered on 
each face of the unit cell. In the ferroelectric tetragonal 
phase, the oxygen octahedron is elongated along [0 1] 
and also shifted by 54 pm with respect to the Pb 2+ posi- 
tion, the titanium ion being displaced from the center of 
cell along [0 1] by 17 pm, respectively. As a result, the 
titanium ion occupies an off-center position within the 
shifted oxygen octahedron (cf. table III Bl figures ^ EJ ■ 

Iron is believed to substitute Ti 4+ as Fe 3+ at the B- 
site in the ABO3 perovskite structure. In order to ver- 
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Figure 3: Comparison of the experimentally obtained axial FS parameter b® = 35.28 GHz (solid horizontal line) for the Fe 3+ 
center in PbTiOj- with calculated values for i>§ obtained from the NSM analysis in dependance of the displacement d from the 
Ti 4+ position. The positions of the oxygen vacancies, of the center of the oxygen octahedra, as well as of the DFT-calculated 
shifts of the Fe 3+ ion are indexed by arrows as reference, (a) - substitution of the Fe 3+ ion either at the A-site. (b,c) - 
substitution of the Fe 3+ ion at the B-site with oxygen vacancies either at Ol or Ol'. The model for 'free iron', i.e. without an 
oxygen vacancy in the nearest neighbor shell, is represented in all models by a dashed line and the model for the iron-oxygen 
vacancy associate by a dash-dotted line in case of static oxygen position and by a solid line for the model using relaxed oxygen 
positions. 



ify this hypothesis, in the calculation we also considered 
the possibility that it substitutes at the A-site for Pb 2+ . 
Furthermore, to decide whether the oxygen vacancy is 
situated in the first coordination sphere or if it is lo- 
cated more distantly, two more structure models were 
investigated, with and without oxygen vacancy in nearest 
neighbor positions. The model of a directly coordinated 
oxygen vacancy can be realized in two different arrange- 
ments, in which the vacancy may substitute either for 
the apical Ol or 01' position (cf. figure QJ. A vacancy 
coordination at the equatorial 02 and 03 positions can 
be excluded as the orientation of the Fe^ — V" defect 
dipole has to be along the c-axis, due to the EPR results 
that indicate axial symmetry at the iron center—. For ori- 
entations of the oxygen vacancy along the a and b axes, 
the FS tensor would be of only orthorhombic symmetry. 
For a final model, the effect of a likely relaxation of oxy- 
gen ion positions in case of a directly coordinated oxygen 
vacancy has been taken into account. The relaxed ionic 
positions were obtained from the DFT calculations (cf. 
section ITTTCj) . 

The position of the oxygen octaeder is defined using 
the crystal coordinate system given in table HH1 Here, 
the oxygen vacancy is located either at (0.5000 0.5000 
0.1239), or at (0.5000 0.5000 1.1239) in fractional coordi- 
nates, and the reference position of iron is at the titanium 



position in the undoped system at (0.5000 0.5000 0.5415). 
For the calculation of b®, the position of Fe 3+ was varied 
along [0 1], the parameter d defining the shift of the 
ion. A positive sign of d defines a shift towards the Ol' 
oxygen, i.e. d = z x 416.8 pm. 

The results are presented in figure |3| depicting the de- 
pendence of calculated axial FS parameters b\ on d for 
the different structural models. In this figure, the ex- 
perimentally obtained value for b® is represented by a 
solid horizontal line. As reference, the positions of the 
corresponding oxygen vacancies Vq', Vq',, the position 
of the equatorial plane of oxygens through 02, 03, as 
well as the positions for the iron as obtained by the DFT 
calculations (cf. section ITlI C|) are indicated by arrows. 

Four conclusions immediately emerge from the NSM 
calculations: as compared to earlier studies for the iron 
center in lead titanat o 9 ' 10 , in which static oxygen posi- 
tions were used, the inclusion of relaxed positions for the 
equatorial oxygens towards the vacancy yields substan- 
tially refined and quite different curves for b®- Further- 
more, the extremum of the b® dependence is no a-priori 
viable position for the iron displacement as has earlier 
been exploited as first gues o 9 ' 10 ; the intersection points 
with the experimental value, not necessarily coinciding 
with the extremum, should be used instead. Finally, 
agreement with the experimental valued can be obtained 
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only with a model, in which Fe 3+ is substituted at the 
B-site with a directly coordinated oxygen vacancy, thus 
forming an Fe^ — V" defect associate. 

In the pure, unsubstituted ferroelectric phase, Ti 4+ is 
displaced considerably from the center of the oxygen oc- 
tahedron (cf. table III Bfl . In contrast, there are two pre- 
dicted positions for the Fei^j" ion with an adjacent oxygen 
vacancy for each of the two models: for the Fe^ ; — Vq" 
associate the iron ion is either displaced by 14 ± 5 pm 
from the Ti 4+ position away from the vacancy towards 
the center of the truncated octahedron, or it is shifted 
considerably by 23 ± 5 pm towards the vacancy. For the 
alternative Fe^ — Vq*, associate, the Fei^j~ ion both times 
is shifted towards the vacancy, where one position is close 
to the plane of equatorial oxygens at 13 ± 5 pm and the 
other position is at 76 ± 5 pm. In order to distinguish 
between these positions, further information as provided 
by the DFT calculation is needed. The finally preferred 
structure (vide infra) of the Fe^ ; — V" defect associate 
from the NSM calculations is presented in figure 0] 



C. DFT calculations 

1. DFT investigations of undoped PbTiOz 

First, the iron-free host lattice of PbTiC>3 was studied. 
Using the VASP code, lattice parameters a — 386.6 pm 
and c — 403.9 pm were obtained, which is in agreement 
with earlier calculations^. These values are smaller than 
the experimental values as expected from LDA calcula- 
tions. The percentage difference between the experimen- 
tal and calculated values are different for the a and c 
direction resulting in a calculated c/a strain of 1.04 com- 
pared to the experimental value of 1.07. The calculated 
Ti-0 bonds are given in figure 0^c) . 

Next, one oxygen vacancy was inserted in the 2x2x2 
supercell without providing charge compensation, i.e., 
the electronic structure of PbgTi 8 023 was calculated for a 
determination of the local relaxation of the ions around 
the vacancy. As shown in figure 0] (d), one finds that 
the nearest neighbor Ti ions move away from the va- 
cancy site, which was also predicted in a previous LDA 
calculation^. If the vacancy is located at the 01 site 
(cf. figure^), the Ti moves towards 01' and the Ti-01' 
distance decreases in the LDA simulation from 224.5 pm 
(in the perfect lattice) to 193.2 pm. The Ti ion clos- 
est to the vacancy in positive c-direction moves from its 
original position by d = +25.4 pm while the nearest Ti 
ion on the opposite side of the vacancy is displaced by 
d = —3.9 pm. If the vacancy is located at the 02 site, Ti 
moves towards the 02' site and the Ti-02' bond distance 
decreases from 195.3 pm to 187.6 pm. 



2. DFT investigations of iron-doped PbTiOz 

The theoretical studies of the system Fe':PbTi03 have 
been performed using an 2 x 2 x 2 supercell, in which 
now Ti is partly replaced by Fe. We have studied differ- 
ent possibilities of ion arrangements for the defect struc- 
tures. In a model series A, we have replaced two Ti by 
Fe', Fe^j, and for charge compensation we have created 
one O vacancy, Vq', that is the supercell has the com- 
position PbgTi6Fe2023- In two other model series, we 
have replaced only one Ti by Fe'. In the model series 
B, we consider the neutral defect [Pb8Ti 7 Fe023] and in 
the model series C, the charged defect [PbsTi 7 Fe023] + , 
charge compensation obtained by assuming a homoge- 
neously charged background. 

The ion arrangement in the three model series A to C 
are as follows. In model Al, the two Fe^ are next nearest 
neighbors and the V" is located between these two Fe^j 
forming a clustered defect oriented in the c direction. 
Model A2 differs from model Al only by the assumption 
that one of the Fe^; is an isolated point defect, whereas 
for model A3, we have assumed that all point defects are 
isolated. Model A4 to A6 vary from model Al to A3 by 
the position of the oxygen vacancy, which is located in 
the a (or b) direction (i.e., at 02 or 03 in figure^) rather 
than in the c direction (01 in figure For the systems 
of model series B [Pb 8 Ti 7 Fe0 23 ] and C [Pb 8 Ti 7 Fe0 2 3] + 
with only one Fexi in the 2x2x2 supercell we have 
again studied the arrangement of V" — Fe^i neighbored 
in a or c direction (models B4 or Bl, B2 and C4 or CI, 
C2) as well as the arrangement with both Fe^; and Vq' 
as isolated point defects (models B3 and C3). For the 
defect associate Vq' — Fe' Ti neighbored in the c direction 
we distinguished between the case where the vacancy site 
is at 01 (models Bl and CI) or at 01' (models B2 and 
C2) as the distances of the 01 and 01' to Fe^j are not 
the same. 

For all these models we have performed a structure 
optimization varying the atomic positions but taking the 
lattice constants a and c fixed to the values obtained 
from the DFT results for the pure host lattice given in 
table III 51 Considering the total energy of the optimized 
structures it is concluded for all three model series A to C 
that the oxygen vacancy in c direction is always favored 
in energy, as compared to an orientation along a. For 
models Al to A6 the clustered defect Fe^; — V" — Fe^i 
along the c direction (model Al) gives the lowest energy 
.Eai followed by model A2 (Fe^j — V" placed along the 
c direction) with energy Ep^.-, where Eki — -Eai = 0.59 
eV per supercell. 

Also for the model series with only one Fe^j defect the 
direction of the Fe^ — V" associate along c (models Bl 
and CI) is preferred over the a direction. More precisely, 
the most stable arrangement is the arrangement were 01 
is missing and not 01'. The models Bl and CI are lower 
in energy than B2 and C2, namely Eb2 — Ebi = 0.35 eV 
and Ec2 - E C i = 0.34 eV. The arrangements B4 and C4 
with an orientation of the Fe^i — Vq* dipole along the 
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Figure 4: Comparison between the different structure models for Fe:PbTi03 deduced from the theoretical and experimental 
investigations. The given numbers indicate atomic distances (in pm) calculated by DFT (cf. section 3.3.2) and the bracketed 
values are experimentally and semi-emprically obtained values. The relative energies for comparable models are shown below 
the figures. 

(a) undoped reference structure obtained from DFT calculations and from synchrotron powder difraction (cf. table ITTTj^l : 

(b) Fexi-V"- Fe'xi complex in Pb8Ti(;Fe2023 oriented along c-axis (model Al); 

(c) Fexi-V" defect associate in [PbsTi7Fe023] + oriented along c-axis where the Ol oxygen is missing (shorter bond, model 
CI). The bracketed values are obtained with Newman superposition model on the basis of measured FS data. This model is 
the energetically most favored among models with one oxygen vacancy. 

(d) Fexi and V" as isolated defects in [Pb 8 Ti 7 Fe023] + (model C3); 

(e) Fexi-V" defect associate in [PbgTi7Fe023] + oriented along c-axis where the oxygen Ol' is missing (longer bond, model 
C2); 

(f) Fe^i — V" defect associate in [Pb8TiyFe2023] + oriented along a-axis where the oxygen 02 is missing (model C4); 
Models B1-B4 which consider defects in charged PbgTi7Fe2023 are in the same energetical sequence. 



a axis is less stable than the favored arrangement within 
the model series by E-q^—E-qi = 0.45 eV and Ec4—Eci — 
0.38 eV. The associated defect Fe^i — V" is favourable 
compared to the isolated defects by E&3 — Em = 0.27 eV 
and Ec3 — Eqi = 0.13 eV. Applying the spin-polarised 
version of the LDA method changes the energy by less 
than 0.02 eV per 2x2x2 supercell and therefore has 
not been pursued any further. 

Next, the movement of the ions around the vacancy 
shall be described. Qualitatively, we find the same trends 
for all models considered here, namely the cation nearest 
to the vacancy moves away from the vacancy site, which 



results in a shortening of the metal oxygen bond opposite 
to the vacancy. 

Comparing the undoped structure (figure 21 (a)) and 
the defect structures (figure 01 (b)-(f)) one sees that the 
position of the metal ion above the vacancy moves up- 
wards above the plane of the four oxygen ions which has 
already been found for Tixi in a previous calculation by 
Park and Chadi—. The displacement of the Fe ion from 
the original Ti position in the undoped host lattice repre- 
sented by the distance d is equal to +46.6 pm and —11.6 
pm for the two Fe' Ti in model Al. For the structures 
with only one Fe^; we get the following values for d, 
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c?A2 = +32.9 pm g?bi = +31.2 pm, dci = +33.5 and 
d c'2 — —6.8 pm. 

Considering the atomic arrangement and relaxation we 
finally point out that even for the pure host the weak O 
- Ti bond of 243 pm compared to the other O - Ti bonds 
of 174 and 198 pm are not described quite satisfacto- 
rily by the LDA results, see figure 0] This also holds 
for GGA (Generalized Gradient Approximation) calcu- 
lations which we have additionally performed using the 
potential of Perdew and Wang 4 ^. Therefore, we have re- 
peated all calculations under the constrain of leaving the 
strain constant, c/a — (c/a) exp — 1.0721 by both LDA 
and GGA. From these additional calculations we get the 
same main conclusions described above, namely that the 



(a) 



associate Fe^ 



V " in 



direction has the lowest energy 



and compared to the original Ti and O positions the Fe^; 
is moving away from the vacancy V". 

The resulting electron states are presented in figure 
in which the density of states (DOS) and the partial DOS 
based on the Mulliken population analysis is plotted for 
PbTi03 (figure (a)). The lowest states in figure 0(a) 
are the 6s-like states of Pb. Above these states we find 
the oxygen 2p-like states mixing with the 6s- and 6p-like 
states of Pb and the 4s/4p-like states of Ti. This indi- 
cates that the bonding has a distinct covalent character 
as found also in previous calculationo 15 i 18 i 46 . The lowest 
states of the conduction bands are predominantly 3c?-like 
states of Ti for pure PbTi0 3 . For Pb 8 Ti 6 Fe 2 023 in figure 
[21(b), additional Is and 2p basis functions centered at the 
vacancy site have been used. No new states are found in 
the band gap, which otherwise would indicate the exis- 
tence of color centers at the oxygen vacancy. Compared 
to figure (a) we get additional Fe like states and one 
can see that the highest occupied states are now the 3d- 
like states of Fe. The oxygen, for which the partial DOS 
are displayed in figure GJb), is neighbored to Fe in the c- 
direction and labelled as O(Fe'). Its 2p z -like states show 
a common peak with the 3d-like states of Fe in an energy 
range of —7 to —6 eV indicating an Fe - O bond. 



IV. DISCUSSION 

The most striking prediction from the NSM calcula- 
tions is that the Fe 3+ functional center in lead titanate 
forms a charged Fe^; — V" defect associate with an oxy- 
gen vacancy in the local oxygen octahedron, the iron be- 
ing substituted at the B-site of the perovskite AB0 3 lat- 
tice. The similarity of ionic radii between Ti 4+ at 68 
pm and Fe 3+ at 64 pm as compared to Pb 2+ with a 
radius of 124 pm supports this assignment. In fully or 
partly ionic compounds, vacancies are charge balanced 
by other defects forming an overall neutral system. It 
can be assumed that partial charge compensation takes 
place at a nearest-neighbor oxygen site in the octahe- 
dron, because the resulting Coulomb interaction is the 
most important driving force for association. This as- 
signment is in accordance with first-principle calcula- 



DOS 

? 


Pb' 

S 


Pb ' 

) ' 


Ti 

S 


Ti 

p 


Ti ' 

> 


01' 

) 


02' 
) 


> 


j 

> 


1 
) 






1 
\ 


i 


) 
1 


10,00 


2.00 


2.00 


2.00 

N(E) 


2.00 

leV 1 ] 


2.00 


2.00 


2,00 




0,F= 



0(F« 



Fe 



Fe 



Fe 



2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 2.00 

N(E) leV] 



Figure 5: Density of states (DOS) for undoped PbTiOs (a) 
and Fe:PbTi0 3 (b). Total DOS and partial DOS of Pb, 
Ti(s,p, d), 01-3(s,p). The horizontal dotted line corresponds 
to the top of the valence band. 



tions that predict the binding of oxygen vacancies to 
acceptor-type impurities 1 ^, as well as to our DFT inves- 
tigations that favor defect associates. Experimental find- 
ings indicate, however, that for other acceptor centers in 
similar compounds, like chromium-doped PbTiO^ and 
copper-modified Pb[Zr .54Ti .46]O3 48 , no such associa- 
tion is present. Alternatively, like in iron-doped SrTiOs, 
an equilibrium between 'free' Fe^j centers and Fe^ ; — Vq" 
defect associates can be present. With respect to the 
ionic mobility of free oxygen vacancies, the iron - oxygen 
vacancy defect dipole complex will be rather immobile in 
the ceramic. Hence, charge transport will be consider- 
ably hindered. Furthermore, the charged defect agglom- 
erate may influence the poling properties by providing 
pinning centers for domain walls^S. Moreover, dipolar 
defect complexes have been demonstrated to be able to 
also pin the polarization of the surrounding crystal 1 ^. 
In particular, oxygen- vacancy related defect dipoles have 
been shown to be involved in voltage offsets leading to 
imprint failure^ and are suggested to play a crucial role 
in electrical fatiguo 19 i 52 i 53 . 

With respect to the reliability of the predictions ob- 
tained by the semi-empirical NSM approach, the follow- 
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ing points should be mentioned. First, crystal distortions 
near the substituted ion and contributions from ions or 
vacancies more distant than neighboring ligands are ne- 
glected. Second, the main assumption of the NSM is that 
the spin-Hamiltonian parameters result from individual 
crystal field contributions of every nearest neighbor ion. 
Finally, the model is based upon the calculation of sin- 
gle electron-derived charge densities, whereas the FS in- 
teraction is related to two-electron expectation values. 
Apparently, the NSM has proven to yield reliable results 
for determining the second- and fourth-rank finestructure 
parameters in S'-state ions, probably because the intrinsic 
parameters have been refined over several decades using 
a large set of experimental data. If iron at a B-site of the 
perovskite ABO3 lattice can be considered as part of this 
set, the structural data can be considered as reliable. 

Within these limits, the results can be interpreted as 
follows: the Fc 3+ ion with the nearby oxygen vacancy 
probably has a position near the center of the truncated 
oxygen octahedron, i.e. at the (0.5000 0.5000 0.6288) 
position, whereas the second position in accordance with 
NSM prediction can be discarded because of DFT-based 
results. 

As mentioned above, for all DFT calculations a 2 x 2 x 2 
unit cell was used. This size of the unit cell seems to be 
too small to avoid interaction of the defects, and the re- 
laxation at the defect site could be smaller if larger unit 
cells are used. However, for computation time reasons 
we could not expand our unit cell. Nevertheless, it is 
most encouraging to find that the local structure with re- 
spect to Fe 3+ positioning, predicted by DFT calculation 
and by the semi-empirical model are quite similar. The 
agreement is surprisingly good and might be fortuitous, 
considering the rather small supercell accomplished by 
an effective 12.5 % iron doping 

The quantum mechanical simulations show, that the 
association of the defects is preferred over the arrange- 
ment as isolated point defects. Additionally, the arrange- 
ments for an orientation of the defect associate along 
the c-axis are more stable than the ones along the a 
and &-axes, confirming the conclusion drawn from the 
axial symmetry of the FS tensor—. The off-center posi- 
tion of the Fe ion in the oxygen octahedron leads to two 
Fe^i — V" distances in c direction and therefore to two 
different positions of the oxygen vacancy. According to 
the DFT investigations the vacancy at the 01 position in 
figure n is energetically favored, in which the long 01'- 
Fe^i distance decreases significantly. 

Finally, because the observed EPR spectra™ can al- 
most exclusively be interpreted in terms of Fe^j — V" 
associates, one can conclude that no 'free' Fe 3+ signals, 
i.e., iron ions without associated V" and hence with 
considerably smaller FS values, are present. Consider- 
ing the condition for overall charge compensation there 
is a charge mismatch, because the Fe 3+ ion substitut- 
ing for the Ti 4+ ions are singly negative charged (Fe^ ; ), 
whereas the associated oxygen vacancy is doubly positive 
charged {V") with respect to the neutral lattice. Hence, 



additional mechanisms for charge compensation have to 
be discussed. Candidates for charge compensation are 
either free electrons (e') trapped in the lattice, lead va- 
cancies (Vp h ), and the formation of positively charged 
cations, such as Pb 2 + -> Pb+, Ti 4 + -> Ti 3+ , or Fe 3+ 
— > Fe 2+ . However, there no evidence for color centers in 
the DFT calculations and the variable valency ions Pb + , 
Ti 3+ are paramagnetic with g-values at about g w 2.0 
to 1.9, which were not observed in the EPR spectra 
and thus can be excluded. Since intrinsic double neg- 
atively charged lead vacancies Vp h have been suggested 
as additional charge compensation^, the overall electro- 
neutrality condition for iron-modified lead titanate thus 
is proposed to be given by 

[Vo} = lFe' Ti } + ±[V P ' h } (7) 

This model is supported by the inherent loss of PbO 
during processing, for which reason a natural intrinsic 
V" — Vp' b di-vacancy pair was proposed^, its existence, 
however, being currently controversially discussed on the 
basis of DFT calculations 2 ^ 3 .. 

In principle, charge-compensation could also be ob- 
tained by the creation of Fe' Ti — Vq*— Fe^; defect asso- 
ciates without the need for lead vacancies at all. The 
existence of such a defect structure is also supported by 
the DFT results, predicting this arrangement as of lowest 
energy. In this structure, the DFT calculations predicts 
two different iron sites, for which reason two Fe 3+ EPR 
spectra with different FS parameters are expected. Using 
the calculated displacements of d = +46.6 pm and —11.6 
pm, the corresponding NSM estimates are 6° = 9-6 GHz 
and 6° = —34.7 GHz, respectively. These values differ 
considerably from the experimentally obtained one, and 
for this reason this model can be discarded. Furthermore, 
for samples with low iron doping, this prediction of an 
Fe^i — V"—Fe' Ti defect associate might not be upheld. 
Finally, from the present EPR experiment there is no fur- 
ther evidence for this hypothesis, because no indication 
for strongly dipolar coupled iron centers is detected in 
the EPR spectra. 

V. CONCLUSION 

In summary, a charged Fe^i — V" defect associate 
in lead titanate has been identified and its microscopic 
structure has been determined based on a comparison of 
FS data with results of semi-empirical NSM and DFT 
calculations. The refined structure comprises informa- 
tion about the structural relaxation around the iron func- 
tional center in lead titanate. In the model presented, 
Fe 3+ is substituted as an acceptor center at the per- 
ovskite B-site with a directly coordinated oxygen va- 
cancy. The position of the iron ion in the ferroelectric 
phase is found to be almost centered in the unit cell, dif- 
ferent from the bulk B-site Ti 4+ ions in the ferroelectric 
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phase that are considerably displaced, but similar to the 
B-sitc positions in the paraelectric phase. 

The orientation of the Fe^; — V" defect dipole is found 
to be along the crystallographic c-axis. This can be con- 
cluded from EPR results which indicate axial symmetry 
at the iron center. Any other orientation of the defect 
dipole would result in an FS tensor of lower than ax- 
ial symmetry. Furthermore, the DFT calculations con- 
firm this assignment by showing that the total energy 
of the arrangement along the c-axis is more stable than 
the ones along the a and 6-axes. Using high-resolution 
synchrotron powder diffraction, a quite large c/a-ratio of 
1.0721 was measured, again pointing to the uniqueness 
of structural relaxation along c. 

The impact of iron-doping in lead titanate on the 
macroscopic piezoelectric properties can now be rational- 
ized from a microscopic point of view, because the iron - 
oxygen vacancy defect dipole complex will be rather im- 
mobile in the ceramic as compared to the ionic mobility 



of free oxygen vacancies. Thus, charge transport will be 
considerably hindered and the charged defect agglomer- 
ate may furthermore influence domain-wall motion and 
poling properties by providing pinning domain walls. 
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